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A theory of spin-polarized transport in inhomogeneous magnetic semiconductors is developed and 
applied to magnetic/nonmagnetic p-n junctions. Several phenomena with possible spintronic appli- 
cations are predicted, including spinvoltaic effect, spin valve effect, and giant magnetoresistance. It 
is demonstrated that only nonequilibrium spin can be injected across the space-charge region of a 
p-n junction, so that there is no spin injection (or extraction) at low bias. 
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Semiconductor spintronics is a developing field where 
active control of spin dynamics is projected to lead to de- 
vice applications with more functionality and better inte- 
grability with traditional semiconductor technology than 
metal-based spintronics fl. Spin injection into a semi- 
conductor has been demonstrated H[| , but an important 
question, what this injected spin is useful for, remains 
largely unanswered (though some semiconductor spin- 
tronic device schemes have been proposed [QH])- This 
Letter addresses the problem of spin and charge bipolar 
(electron and hole) transport in inhomogeneously doped 
magnetic semiconductors (such as BeMnZnSe, CdMnSe, 
or GaMnAs, whose carrier g-factors are large due to 
magnetic doping). We show, as an important conse- 
quence of our theory, that a magnetic/nonmagnetic p-n 
junction has significant potential for device applications, 
and predict large magnetoresistance for such a junction. 
We also show that the magnetic junctions display spin- 
voltaic and spin valve effects: a current, whose direction 
changes with the direction of applied magnetic field or 
injected spin polarization, can flow without any applied 
bias. These effects should be useful for sensing magnetic 
fields, and for probing spin polarization and spin relax- 
ation. 

To introduce the equations for spin and charge bipo- 
lar transport in inhomogeneously doped magnetic semi- 
conductors, consider a semiconductor doped with iV a (r) 
acceptors and Nd{r) donors, and with magnetic impuri- 
ties, whose density varies in space and whose presence 
leads to large <?-factors for electrons and holes, g n {r) and 
g p (r). In a homogeneous magnetic field B the carrier en- 
ergies are Zeeman-split: spin up (A = 1 or f ) and spin 
down (A = —1 or J.) electrons have their energy shifted 
by — Xq( n (r) = — A<7„(r)//B-B/2, where \xb is Bohr mag- 
neton and q is the proton charge; the energy of holes 
changes by Xq( p (r) = Xg p (r)figB/2. Carrier charge cur- 
rent densities, resulting from the electric field E = — V0, 
nonuniform magnetic "potentials" £, electron (hole) den- 
sities n (p), are 

JnA = qHn\n\E + qD nX Vnx - gA^„An A VC„, (1) 
J p a = qn P \P\E - qDpxVpx - q\pp\p\VC,p, (2) 

where u and D stand for mobility and diffusivity, and 



the "magnetic" drift q\7( forces carriers with opposite 
spins to go in opposite directions. The above equations 
follow from the requirement that carrier densities in an 
inhomogeneous environment have the quasiequilibrium 
form, for example, n\ = nx(q<fi+Xq(, n +ftx), where n\{fi\) 
is the number of electrons per spin A at (nonequilibrium) 
chemical potential fix , and the expression for the current 
Jn\ = n\p\VjI\ M. From Eqs. [I] and ||, the carrier 
charge and spin current densities J = Jj + Jj and J s — 
J-f — J i are: 



J„ = a n E - <7 sn VCn + qD n Vn + qD sn Vs 7 , 
J p = (T P E - cr sp VCp - qDpVp - qD sp Vsp, 



and 



J s « = c s „E - c„VC„ + qD sn \7n + qD n Vs n , 
J sp = c sp E - CpVCp - qD S pVp - qDpVsp. 



(3) 
(4) 

(5) 
(6) 



Here n = n-f + n± and s„ = n-j- — n^, and we in- 
troduced electron charge and spin conductivities <j n = 
q(fi n n + u s „s„) and a sn = q{p sn n + (J. n s n ), where //„ = 
(AM + AM)/2 and fx sn = (^i nT - am)A anci similarly 
for diffusivities. Analogous notation is used for holes. 
Equations ^-|[ which are a generalization of the Johnson- 
Silsbee magnetotransport equations 0, reflect the spin- 
charge coupling in bipolar transport in inhomogeneous 
magnetic semiconductors: a spatial variation in spin den- 
sity, as well as in gV£, can cause charge currents, and 
spin currents can flow as a result of a spatial variation of 
carrier densities and <f>. 

Generation and recombination of electrons and holes 
are assumed to be mostly due to band-to-band processes. 
Furthermore, electrons (holes) with a given spin are as- 
sumed to recombine with holes (electrons) of either spin. 
The stationary continuity equations for electrons and 
holes read 

As r 



V7 J » A I / \ | ^X 
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(7) 
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Here no and po are the equilibrium carrier densities, w 
is the (generally spin dependent) band-to-band recombi- 
nation rate, and T\ is the spin relaxation time. Spin 
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FIG. 1. Band-energy schemes for magnetic/nonmagnetic 
p-n junctions with magnetically active electrons (arrows). 
Holes (circles) are unpolarized. Electrons from the magneti- 
cally active n-region (A) (discernible by the conduction band 
split) can be injected into the nonmagnetic p-region only 
at large bias; similarly for spin extraction from a magnetic 
p-region (B), where electrons are minority carriers. If spin 
of the majority electrons is out of equilibrium (they are spin 
polarized), and the p-region is magnetic (C), a giant magne- 
toresistance and spinvoltaic effects arise. The former can be 
observed in a scheme (D), where electron spin is injected from 
a magnetic heterostructure N into the nonmagnetic n-region, 
which forms a p-n junction with a magnetic p-region. 

relaxation equilibrates electron spins while preserving 
nonequilibrium electron density, so, for a nondegener- 
ate semiconductor, s n = a n on, where a„o = s n o/no = 
tanh(C„/Vr) is the equilibrium electron spin polarization 
(Vt = ksT/q, with ks being the Boltzmann constant 
and T temperature). In general, s n = n[p-\- qCn) ^^{(1 — 
qCn), where p is obtained from the known nonequilibrium 
electron density n = n(jl + q^ n ) +n(p, — qC, n )\ similarly for 
holes. To include transient effects, terms —dn\/dt and 
dpx/dt need to be added to the left sides of Eqs. and 

Equations |l| and ^ (or, f|j(|), and |[ and Poisson's 
equation V • eE = p, where p = q(p — n + Nd — N a ) and 
e is the dielectric constant of the semiconductor, deter- 
mine the distributions of charge and spin in a magnetic 
semiconductor under applied bias V. For ferromagnetic 
semiconductors ( needs to be evaluated self-consistently 
(similarly to </>). Indeed, the magnetization depends on 
the carrier density (as the mobile carriers mediate the 
magnetic interactions), and the carrier density depends 
on magnetization, as magnetization determines £. 

We have made a number of simplifying assumptions 
in our theory which are not strictly valid in real mag- 
netic semiconductors such as GaMnAs. For example, 
electron and hole band states are treated as simple spin 
doublets, although the band structure is usually more 
complicated in the presence of spin-orbit coupling. We 
also do not take into account the fact that magnetic im- 
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FIG. 2. Calculated electron (n), hole (p), and spin (s) den- 
sity profiles for a magnetic p-n junction with q£ = 0.5fcsT 
in the (magnetic) n-region, under forward applied bias of 
V — 0.8 Volt. The thin dashed line indicates the doping 
profile Nd — N a (not to scale). 

purities change the gap and can offset the bands, even 
at B=0. Finally, in some magnetic semiconductors (even 
nonferromagnetic) carrier densities may strongly affect 
the response of the magnetic ions to magnetic field, in 
which case a selfconsistent calculation is needed. In prin- 
ciple, our formalism can be tailored to include these more 
complicated effects, but this can only be done on a case- 
by-case basis (if the effects are well understood, which is 
currently not the case). Our approach here is to con- 
sider a simple but reasonable model to illustrate new 
phenomena that can occur in inhomogeneous magnetic 
semiconductors . 

We now apply the above theory to the problem of a 
magnetic/nonmagnetic p-n junction, with only electrons 
being magnetically active (£ p = 0) ||. Two cases are 
considered (Fig. |l|) : (i) the magnetic n-side to study spin 
injection, and (ii) the magnetic p-side, to study spin ex- 
traction and magnetoresistance phenomena. We assume 
perfect ohmic contacts (both carriers and spins at equi- 
librium) as boundary conditions, except when we study 
the dependence of the current on nonequilibrium spin (in 
which case we keep the ohmic contacts for carrier densi- 
ties only, not for the spin). We only deal with nondegen- 
erate semiconductors, as these are analytically tractable, 
and yield the largest magnetoresistance. 

To be specific, we consider a L = 12 /im long GaAs p-n 
junction at room temperature, doped with N a = 3 x 10 15 
cm' 3 acceptors to the left, Nd — 5 x 10 15 cm~ 3 donors 
to the right j5|, and with magnetic impurities inducing 
electronic (/-factor which follows the shape of N a (Nd) for 
the magnetic p (n) region. Holes are unpolarized, so we 
omit the label n when dealing with spin-related quan- 
tities, and reserve that symbol for the n-region-related 
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FIG. 3. Calculated spin polarization profiles for different 
forward bias and spin relaxation rate. Spin injection (top) 
from the magnetic n-side into the nonmagnetic p-side occurs 
only at large bias (indicated by the numbers in Volts). The 
largest injection in the graph is for V = 1.5 Volt, with spin 
relaxation time of lOOTi . Similar behavior is observed for spin 
extraction from the nonmagnetic n-region into the magnetic 
p-region (bottom). The magnetic splitting is g£ = 0.5ksT. 

variables. Figure [l] depicts our geometry and the distri- 
bution of electrons and holes for various cases considered 
below. The materials parameters || are the electron and 
hole diffusivities, D n = 10D p = 103.6 cm 2 -s _1 , and mo- 
bilities fj, n = lOjUp = 4000 cm 2 'V -1 -s -1 . The intrinsic 
(nonmagnetic) carrier density is rii — 1.8 x 10 6 cm~ 3 , 
the dielectric constant e = 13.1 of the vacuum permit- 
tivity; the calculated built-in- voltage is 1.1 V at B = 0. 
Recombination rate w — (1/3) x 10~ 5 cm 3 -s _1 , and the 
spin relaxation time T% = 0.2 ns. The minority |^Cj[| dif- 
fusion lengths are Q L n — 1 /im, L. p = 0.25 /im, and 
the electron spin diffusion length in the n (p) region is 
L sn = 1.4 /im (i sp =0.8 /im). Figure || illustrates the 
doping profile and a typical profile of carrier and spin 
densities. 

We first ask the important question whether spin can 
be injected and extracted into/from the nonmagnetic re- 
gion. Figure || shows the results of our numerical calcu- 
lations, where we plot spin polarization a = s/n (note 
that we are interested here in spin polarization of den- 
sity, not current ||). At small bias (below the built-in- 
value), which is also the small injection limit, there is 
no significant spin injection or extraction. As the bias 
increases, the injection and extraction become large and 
intensify with increasing T±. The reason why there is no 
spin injection (and, similarly, extraction) at small bias 
is that although there are exponentially more, say, spin 
up than spin down electrons in the magnetic side, the 
barrier for crossing the space-charge region is exponen- 
tially larger for spin up than for spin down electrons (see 



Fig. |l|). Those two exponential effects cancel each other, 
so there is no net spin current flowing through the space- 
charge region. The numerical aproach is indispensible for 
obtaining the high-injection results, where we observe a 
large nonequilibrium spin polarization around the space- 
charge region, resulting in spin injection and extraction. 

The current through a magnetic/nonmagnetic p-n 
junction depends on magnetic field. This dependence 
has two sources. First, a magnetic field makes the band 
gap spin dependent, leading to the increased density of 
the minority carriers. For example, the equilibrium den- 
sity of electrons in a magnetic p-region (Fig. |l|C) is uq = 
nf cosh((/Vr)/N a (in equilibrium, nopo = n \ cosh(£/Vr) 
from Boltzmann statistics). The current through a p-n 
junction is proportional to the density of the minority 
(not majority) carriers, and so it will grow exponentially 
with B for £ > Vr- This does not happen for unipo- 
lar transport (a homeogenous sample) , where the current 
depends on the density of the majority carriers, which is 
independent of B, as it is fixed by the number of ion- 
ized impurities. The exponential magnetoresistance for 
the system in Fig. [i]B is shown in Fig. |J. The effect 
diminishes with increasing bias. The second way a cur- 
rent can depend on B is if a nonequilibrium spin with 
polarization a(L) is introduced (optically or by spin in- 
jection) into the sample (Fig. |l|C and D). The exponen- 
tial magnetoresistance becomes giant (meaning, the re- 
sistance changes when the orientation of the magnetic 
moment in the magnetic region with respect to the ori- 
entation of the injected spin in the nonmagnetic region, 
changes) as seen in Fig. ^. The current is exponentially 
sensitive to both B and a(L), being large (small) if they 
are of the same (opposite) sign. 

To explain analytically our numerical results, we gen- 
eralize the well-known Shockley approximation H{ for 
the p-n junction to include nonequilibrium spin, which is 
given, in the quasi-neutral regions, by the spin-diffusion 
equation 

D n V 2 s = w{s P - s Po) + (s - ~s)/T x , (9) 

where s = aon, of which n is the solution of the cor- 
responding carrier diffusion equation. We consider the 
magnetic p-region, and impose a nonequilibrium spin po- 
larization a(L) at the right boundary (see Fig. [j]C, where 
also the boundary points used below are indicated). We 
postulate as the second boundary condition for spin, in 
the spirit of the quasiequilibrium approximation j(| , that 
spin current vanishes at x n , the n-side boundary of the 
space-charge region. The spin polarization in the n- 
region is 

a(x) — a(L) [cosh(n) — tanh(n„) sinh(n)] , (10) 

where n = {L - x)/L sn and n„ = {L - x n )/L sn . At 
x = x n , a n = a(x n ) = a(L)/cosh[(L - x n )/L sn ]; it 
vanishes if a(L) — 0, in line with our result of no 
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FIG. 4. Calculated I-B (top) and I-a (bottom) charac- 
teristics for the magnetic p-n junction with the magnetic 
p-region, under different bias (indicated in Volts). We use 
B* = 2fcsT / g\iB (~ 900/g Tesla at room temperature) as the 
scale for magnetic field (B/B* = qC/Vr)- In the top graph 
the red line is for a(L) = 0, and the violet (cyan) line for 
a — 1 and spin relaxation time of 100 (10 4 ) T\. The bottom 
graph is for B/B* = 3, and the red, violet, and cyan lines are 
for the spin relaxation times of T x , IOOTl, and 10 4 Ti. The 
dashed lines are the analytical calculations. In both graphs, 
at V = 1.5 Volt, lines for 10021 and 10 4 Ti almost coincide. 

spin extraction at low bias. In the p-region cv(0) = ocq 
(ohmic boundary condition for spin) and the second 
boundary condition, that at x — x p , the p-side bound- 
ary of the space-charge region, which can be obtained 
from the usual quasiequilibrium condition of constant 
chemical potentials (applied to nondegenerate statistics) 
from 



{olq + a n )/(l + a a n ). The 



nonequilibrium minority carrier density at x = x p is 
n(x p ) = ?io(l + Oioot n ) exp(I//Vr). With these boundary 
conditions the carrier and spin profiles in the magnetic 
p-region are 



n(x) 
n 

so 



= 1 



sinh(a;/iy n ) 
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n(x) smh(x/L sp ) a 
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n 



sinh (xp/L sp ) a 



(1 



,V/V T 



(11) 
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Having the profiles, we calculate the current in the 
magnetic p-n junction. We distinguish equilibrium-spin 
electron J„ and hole J p currents, and noncquilibrium- 
spin-induced current J' n , so that the total charge current 



is J 



J p + J' . The individual contributions are 
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(13) 
(14) 



(15) 



Figure |5| shows how well this analytical model compares 
with our numerical calculation. The small discrepancy 
is caused by the neglect of the recombination processes 
inside the space-charge region; the model breaks down at 
large bias, where the Shockley approximation becomes 
invalid 0. We note that if the spin is injected by the bias 
contact at x — L, one needs to consider a model in which 
a(L) depends on J (since it then obviously vanishes at 
J = 0), and solve Eqs. p~3| ~ |l5| for J. Otherwise the above 
equations describe either the case of an independent spin 
injection, or the region of J where a(L) is independent 
of J. 

The large exponential magnetoresistance effect, J ~ 
exp(|£|/Vr) at large |£| (or B), comes from the increase of 
the minority electron population with |£|. However, once 
a nonequilibrium spin population (finite a n ) is main- 
tained with its sign fixed in the space-charge region, 
a giant magnetoresistance (GMR) should be observed, 
J ~ exp(C/^ T ) at large |C| and a n -> ±1. The GMR 
coefficient is then exp(2£/Vr). In addition to sensing 
B, these effects can be used for an all-electrical probing 
of the injected spin polarization a(L) and spin diffusion 
length L sn , as they both determine a n (a similar device 
scheme is proposed in Fig. [I] D). The key to employing 
the large magnetoresistance effects in spintronics is the 
development of materials with large (/-factors. Even for 
g « 100, the GMR coefficient at T = 100 K and B = IT 
is about exp(0.65), which is close to 200%. A magnetic 
impurity that would give g « 500 would yield 2,500 %, 
(300% at room temperature). The device potential of 
magnetic p-n junctions is enormous, but to take full ad- 
vantage we need materials with large g-factors at room 
temperature [p]. 

Equations [I3|jl5| reveal another interesting phe- 
nomenon, a spinvoltaic effect. While both J n and J p 
vanish for V — 0, J' n stays finite, if a n ^ 0. A current 
flows at zero bias! This is an analogue of the photovoltaic 
effect, where a reverse current can flow by introducing 
photogenerated carriers into the space charge region. In 
the spinvoltaic effect both reverse and forward currents 
can flow, depending on the relative orientation of B and 
a n , so a magnetic p-n junction can act as a spin valve. 
The physics of the spinvoltaic effect is that nonequilib- 
rium spin in the space-charge region disturbs the bal- 
ance between the generation and recombination currents 
(Fig. 0C). If ( > 0, and more spin up electrons are present 
at x n (a n > 0), the barrier for them to cross the region 
is smaller than the barrier for the spin down electrons, 
so more electrons flow from n to p than from p to n, and 
positive charge current results. If there are more spin 
down electrons at x n (a n < 0), the current is reversed. 
We note that, as at V = no spin can be injected by 
the ohmic contact into the n region, the nonequilibrium 
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spin population has to be injected either optically (at a 
distance larger than L n but smaller than L sn from the 
space-charge region), or electronically transverse to the 
flow of J (and also within L sn of the space-charge region). 

In summary, we have explored bipolar magnetotrans- 
port in magnetic/nonmagnetic p-n junctions, which are 
of great potential for emerging semiconductor spintron- 
ics. We have used a simple model to describe the new 
physics that results when the magnetotransport is domi- 
nated by minority, rather than majority charge carriers. 
In particular, we find that spin injection is suppressed at 
low biases, because the spin-dependent barrier to trans- 
port in the depletion layer negates the effect of spin polar- 
ization of the carriers. Furthermore, we predict exponen- 
tial magnetoresistance resulting from the magnetic field 
dependence of the population of the minority carriers, 
and giant magnetoresistance, appearing if the polariza- 
tion of the majority carriers is disturbed externally (by a 
third-terminal spin injection, for example). We also find 
that a nonequilibrium spin-polarized population can lead 
to a charge current-a spinvoltaic effect. Finally, we have 
shown that at small biases the effects of interest can be 
described by a modified Shockley model. 

This work was supported by DARPA and the US ONR. 



[9] Chosen to have typical values for GaAs with no special 
significance about the numbers we use. 

[10] Term minority (majority) referes to the relative popula- 
tion of electrons and holes. 

[11] At low temperature, this is already acheived, g>500 in 
(Cdo.95Mno.o5Se):In. T. Dietl, in Handbook of Semicon- 
ductors Vol. 3, edited by T. S. Moss and S. Mahajan, p. 
1279 (Noth-Holland, New York, 1994). 
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